Critical Dynamics of Dimers: Implications for the Glass 

Transition 

Dibyendu Das 

Indian Institue of Technology Bombay, Powai, India 400076 

Greg Farrell, Jane' Kondev and Bulbul Chakraborty 
Martin Fisher School of Physics, Brandeis University, 
Mailstop 057, Waltham, Massachusetts 02454-9110, USA 

Abstract 

The Adam-Gibbs view of the glass transition relates the relaxation time to the configurational 
entropy, which goes continuously to zero at the so-called Kauzmann temperature. We examine this 
scenario in the context of a dimer model with an entropy vanishing phase transition, and stochastic 
loop dynamics. We propose a coarse-grained master equation for the order parameter dynamics 
which is used to compute the time-dependent autocorrelation function and the associated relaxation 
time. Using a combination of exact results, scaling arguments and numerical diagonalizations of 
the master equation, we find non-exponential relaxation and a Vogel-Fulcher divergence of the 
relaxation time in the vicinity of the phase transition. Since in the dimer model the entropy 
stays finite all the way to the phase transition point, and then jumps discontinuously to zero, we 
demonstrate a clear departure from the Adam-Gibbs scenario. Dimer coverings are the "inherent 
structures" of the canonical frustrated system, the triangular Ising antiferromagnet. Therefore, our 
results provide a new scenario for the glass transition in supercooled liquids in terms of inherent 
structure dynamics. 

PACS numbers: 
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INTRODUCTION 



A variety of systems such as supercooled liquids, colloids, granular matter and foams, 
exhibit a transition from a flowing fluid phase to a frozen solid phase. Jamming due to spatial 
constraints imposed on the elementary constituents of these materials has been proposed as 
a possible common cause of this dynamical arrest ]} 2>EJ- 

Model systems, such as hard spheres, have an important role to play in the investigation 
of such a scenario since they allow for a precise definition of jamming ^ . They are also useful 
in elucidating the precise relationship between thermodynamics and dynamics in materials 
exhibiting a jammed phase p. The entropy-based Adam-Gibbs theory fjrelates the viscosity 
(a dynamical quantity) to the configurational entropy (S' con /)(a thermodynamic quantity) 
through 7] = r]oexp(A/TS con f). The ideal glass transition is associated with the Kauzmann 
temperature at which the configurational entropy vanishes /7]. In this paper, we explore the 
connection between dynamics and thermodynamics in a lattice model of dimers with an 
entropy-vanishing phase transition. 

The dimer model is one of the working horses of statistical mechanics. It provides an 
example of a jammed system which has the added advantage of being exactly solvable 
States of the dimer model are specified by placing dimers on the bonds of the lattice so 
that every lattice site is covered by exactly one dimer; see Fig^ These dimer coverings are 
"locally jammed" ^ as every dimer cannot move to an empty, neighboring bond, without 
violating the packing constraint. Moves that involve loops of dimers and adjacent empty 
bonds, on the other hand, are allowed. An example of such a move for the hexagonal lattice 
involving an elementary plaquette is shown in Fig. Stochastic dynamics of the dimer 
model on the square lattice based on these elementary moves were considered by Henley . 

Most states of the dimer model allow for elementary moves; an example of one which 
does not is shown in Fig. ^p. The smallest move in this case involves a system spanning 
loop, and we call this state "maximally jammed". If we define an energy functional on the 
space of dimer coverings which favors the maximally jammed state, a transition into this 
state can be affected as the temperature is lowered. The central question we address in this 
paper is: What happens to relaxation time scales of the dimer model as the transition to the 
maximally jammed state is approached? We will show that the relaxation is dominated by 
entropy barriers and is sensitive to equilibrium fluctuations near the phase transition point. 
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FIG. 1: (a) Dimer covering of the honeycomb lattice with an elementary loop update indicated by 
the arrows. The numbers are the heights of the equivalent interface, (b) An ordered, maximally 
jammed dimer covering; the equivalent interface is tilted with maximum slope. 



We consider an energy functional that exhibits a continuous transition to the maximally 
jammed state along a metastable line. We find a strong departure from the canonical 
critical-slowing-down scenario [10], which we attribute to the presence of entropy barriers. 
Barriers can be traced directly to the non-local nature of the dynamical moves allowed by the 
jammed states. The longest relaxation time-scale is found to diverge exponentially following 
a Vogel-Fulcher-like form. This is reminiscent of what is observed in fragile glass formers 
Exponential time-scale divergence (activated scaling) is associated with critical points 
in models with quenched disorder 12] and it has been argued that real glasses belong to 
the universality class of random Hamiltonians with such exponential divergences 13| . The 
current model provides an explicit example of a model without quenched disorder which 
exhibits activated scaling. It should be mentioned that the Vogel-Fulcher law has been 
observed in models with entropic barriers [lJ], with traps Q], and within effective medium 
theory [ifj], none of which have an explicit critical point. 



DIMER MODEL 

We consider the dimer model on the 2-d hexagonal lattice of linear size L, having 2L 2 
sites and 3L 2 bonds, with periodic boundary conditions jl^ |. A useful representation of the 
dimer model is given by the height map which associates a discrete interface h(x, y) with 



every dimer covering 



The heights of the interface are denned on the vertices of the dual 



triangular lattice. The height difference A between two nearest neighboring sites is -2 or +1 
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depending on whether the bond of the honeycomb lattice that separates them is occupied 
by a dimer or not; see Fig. Directions in which the height change is +A are specified by 
orienting all the up pointing triangles of the dual lattice clockwise. 



The dimer model has an extensive entropy. The ensemble of equa 



weighted dimer cov- 



erings maps to a rough surface with a gradient- square free energy [18 1. Fluctuations of the 
surface are entropic in origin. A phase transition can be induced in the dimer model by 
including an energy functional which is minimized by a dimer covering corresponding to a 
smooth, maximally tilted surface which corresponds to the maximally jammed state shown 
in Fig. [TJd. 

For periodic boundary conditions the tilt vector, (A x h, A y h), where A XjV h is the average 
height difference in the x or y direction, has only one independent component p 19]. In 
terms of p, the energy functional we consider can be written as: 

W = -^(1 + 8A (1) 

where p, is a dimensionless coupling, proportional to inverse temperature (/3 = 1/kT), that 
drives the transition. 

The entropy of the dimer model as a function of p was calculated exactly H,Q|: 

S( p ) = L 2 { — (l-p) + - f ¥l ' P) dx\n[cosx}\ (2) 
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This function has a maximum at p = which is the e quili brium value at p = 0. For finite 
p this dimer model was previously considered in Ref. [21(. A dimer model with a similar 

HLse transition but with an energy functional linear in p was solved exactly by Kasteleyn 

In the dimer model with the free energy (3F = f3E — S, and the energy and entropy 
given by Eqs. and El there is an interesting phase transition along the metastable line, 
when the order parameter is confined to the free energy well around the zero-tilt state. 
Namely, at //* = it/ (8a/3), the end-point of the metastable line, the order parameter p has a 
discontinuous jump from to 1, characteristic of a first-order transition. At the same time, 
as p* is approached from below, fluctuations of p around diverge, as would be expected 
at a critical point. This transition was discussed in detail in Ref. [21]. In this paper we 
investigate the dynamics of the dimer model near this phase transition point. 
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COARSE-GRAINED DYNAMICS 



As mentioned in the introduction, the hard constraint of no overlapping of dimers, gives 
rise to nonlocal dynamics. We consider stochastic, Monte-Carlo dynamics based on loop 



221. Since 



updates with loops of arbitrary size; a concrete implementation is given in Ref. 
we take periodic boundary conditions, loops with different winding numbers can be formed. 
We restrict loop updates to loops with winding numbers (0,0), (1,0) and (0, 1), only. The 
microscopic transition rates for loop updates are given by Metropolis rules that follow from 
the energy function, Eq. [T] 

Given the microscopic loop dynamics, which satisfy conditions of ergodicity and detailed 
balance, we ask what are the coarse-grained dynamics of the order parameter, p. Since 
the energy function in Eq. ^depend on the global tilt p only, it follows that all updates of 
topologically trivial loops (i.e. those with (0, 0) winding number) have AE = 0. Only when 
system spanning loops with nonzero winding numbers are updated does the energy of the 
state change. This feature naturally leads to fast and slow processes in the Monte-Carlo 
dynamics. On a faster time scale, non-winding loops are updated with no effect on the 
overall tilt of the surface, while on a much slower time scale, winding loops are updated 
causing a change in the tilt of the surface. 

The coarse-grained dynamics of global tilt changes are described by a master equation 
for the probability (P p ), that the dimer model has tilt p, 



dP, 



p 



( j f [_-• p-j./^ . ■■rrx,M J , Hi Pp + Wp !P -i/LPp-l/L + W p> p + i/lP p+ x/l . (3) 

The rates in this master equation obey the detailed balance condition: W p -.i/l >p /W P)P -i/l = 
exp[— (F(p — 1/L) — F(p))}. The usual way of achieving this balance which leads to normal 
diffusive dynamics is to partition the rates symmetrically with W p -\/l, p — exp[— (F(p — 
1/L) - F{p))/2] and W PtP _ 1/L = exp[(F(p - 1/L) - F(p))/2]Q- Equation however, 
features an unusual form for the transition rates between different tilt states. Namely, the 
rates of transitions from higher into lower tilt states (increasing energy transitions) are 
determined by the energy change alone: 

W p ^ 1/L>p = r oe -<*C^-sM) ; (4) 

here To is a constant. This follows from the observation that in order to lower the tilt and 
increase the energy, a system spanning loop, which is always present in a state with p ^ 0, 



needs to be updated. This form of the rates for energy-increasing transitions in conjunction 
with the detailed balance condition implies that the rates of transitions to higher tilt states 
(energy lowering transitions) must be determined by the entropy change: 

w p , p . 1/L = r oe -cs(P-Vi)-5w) . (5) 

The form of the transition rates that we are ar guin g for here, was directly observed in 
numerical simulations of the three coloring model |2J] , which is a close relative of the dimer 
model. The two are equivalent if, in the dimer models, a weight of 2 is attached to each 
loop formed by bonds that are not covered by dimers. 



RELAXATION TIME-SCALES 



The first consequence of the above form of the transition rates is that the time scale 
of relaxation out of a state with tilt p, r p = 1/(W / p _i/l iP + W p+ i/L tP ), is a non- decreasing 
function of p. The exact expressions for r p (measured in units of Tq 1 ), 

-l = e -f Pl iL + e - i [|i n 2+|Mco S (f-Y)]] ; ( 6 ) 

follows from Eqs. 0] and El and it is plotted in Fig. [2^). This time-scale increases monotoni- 
cally with p[25], as in the hierarchical models of Palmer et al. It is in sharp contrast with 
canonical Langevin dynamics around the equilibrium state, for which the time to relax out of 
a macro-state decreases the further the order parameter is away from its equilibrium value. 
(For example, in the Ising model with Glauber dynamics and in the disordered phase, the 
relaxation time out of a given magnetization state decreases with increasing magnetization.) 



AUTOCORRELATION FUNCTION 



To quantify the tilt dynamics we compute the tilt-tilt autocorrelation function C(t), 
defined as: 

C(t) = ~ <"<°»', (7) 
<p(0) 2 } - <p(0)} 2 

with the average taken over different histories of p. An approximate form for the autocor- 
relation function is: 

C(t) e>-(p»v™ ' (8) 
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FIG. 2: (a) The time scales for relaxing out of different tilt states p in the dimer model (scaled 
by L), for a value of \x below the transition, (b) The tilt-tilt autocorrelation function of the dimer 
model. The full line is obtained from Eq. |H] while the dashed line is a result of the saddle point 
evaluation of Eq. |S| Here L = 4096 and time is measured in units of Yq 1 . 
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i.e., C(t) is an equilibrium weighted average of relaxations out of different p states. This ap- 
proximation is based on the assumption that eigenf unctions of the rate matrix are localized 
in p-space and eigenfunctions corresponding to different eigenvalues do not have significant 
overlaps. We will justify this assumption a posteriori by examining the eigenfunctions ob- 
tained from numerical diagonalizations of the rate matrix W PiP >. The asymptotic decay of 
the autocorrelation functions can be extracted by performing a saddle point analysis of the 
sum in Eq. |H] and using a quadratic approximation for the entropy (Eq. |2J. These saddle 
point solution is compared in Fig. Et>) to the result obtained from the sum (Eq. EJ). 
In the limit of p — > p* and t — > oo, saddle point analysis yields: 

C(t) ~ exp{-A(^)[ln(-^-)] 2 }, (9) 
oz p^ p^ p 

showing that C(t) in the dimer model has a log- normal form implying a slower than expo- 
nential decay. From Eq. El we also conclude that the relaxation timescale, r, for the decay of 
C(t) to an arbitrary constant Co, diverges exponentially as p — > p*. This is a Vogel-Fulcher 
type behavior (since p is proportional to /3 — 1/ kT) observed in many fragile glass formers. 
First order corrections to Eq. El lead to an even more rapid increase of time scales, with 
r/lnr diverging as Vogel-Fulcher. 

The coarse grained dynamics defined by the transition matrix elements, Eqs. 0] and El 
were argued to follow from the nonlocal loop dynamics of the dimer models. From this 
form of the ^-matrix all the conclusions about critical dynamics of the dimer model are 
derived. We have confirmed this picture in considerable detail in simulations of the three 
coloring .nodel HQ, which, as discussed easier, is the loop weighted dimer model. The 
loop weights are not expected to affect the qualitative features of the energy and entropy 
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functionals. Indeed, the measured r p for the three-coloring model compare very well 
to the analytical form plotted in Fig. El The numerical evidence for Vogel-Fulcher type 
divergence of the relaxation time scale in this model was reported previously j^. 

The dynamical behavior of the dimer model can be traced back to the interplay between 
the free energy and dynamical barriers. The transition rates presented in Eqs. H]andEl can 
be interpreted in terms of a barrier[2^ B{p) = e (s{p-i/L)~s{ P )+(E{ P -i/L)~E{p)))/2 dividing the 
usual Metropolis rates defined in terms of the free energy: 

W p . l/L , p = rae-PWQ-'W'/Bip) (10) 
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FIG. 3: Barrier height, B(p) (dimensionless) shown as a set of solid lines, and the quadratic 
approximation to the dimensionless free energies of the dimer model (dashed line); p is chosen 
close to /i* and L = 24. Note the logarithmic scale for the barrier height. 

and 

W p , p . 1/L = T e^-^- F ^/B(p) (11) 

The barriers increase exponentially with p as illustrated in Fig. EJ Dynamics of the order 
parameter can be viewed as relaxation in the free energy well in the presence of these barriers. 
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SCALING ANALYSIS OF THE MASTER EQUATION 



The emergence of the Vogel-Fulcher law in the dimer model, based on the master equation 
with transition rates defined in Eqs. El and ^2 f°U° ws from a scaling argument. 

At the critical point, p = p*, the free energy difference between different tilt states close 
to p = vanishes. In this limit, the transition rates are symmetric and given by: 

W p . l/L , p = W p , p _ 1/L = T /B(p). (12) 

The diffusion constant in this symmetric case can be shown to be given by 29 



D(L) = LVo/^^BiPi)) (13) 

where L is the system size and pi = i/L. The longest timescale in the problem is given by 

r(L) = I?/D{L) (14) 

In the limit of large L, the summation in Eq. EH can be replaced by an integral. If we 
make use of the quadratic approximation to the entropy (Eq. |2J) then B(p) = e ( - 8 ^ 3 ^^ +fl *^ Lp 
and the integral can be evaluated analytically, with the result: 

r(L) = A ( e ¥^ _ i) (15) 
1 o 

Note that the same result can be obtained by replacing the summation by the largest barrier 
which occurs at p — 1. The longest time-scale in the system is, therefore, seen to diverge 
exponentially with system size. If all the barriers were equal to one, then we would have 
D(L) = T and t(L) = L 2 /T , which corresponds to simple diffusion. In the presence of the 
barriers, D(L) goes to zero exponentially and this leads to an exponential divergence of the 
relaxation time-scale. 

We can now use a scaling argument to deduce the behavior of the relaxation time-scale 
for p, < /i*. The effective scale (in p-space) over which the free-energy well is flat, and 
therefore the transition rates are symmetric, diverges as the phase transition at p, = //* is 
approached. We argue that this length scale, given by l{p) = 7/^=zp provides a cutoff to 
the summation (or integral) involved in calculating the diffusion constant 

D{1) = WCE^mmBQh)) (16) 
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The longest time scale therefore scales as 



r^) = r{l^)) ] 1(jm)<L (17) 

and 

r(/i) = r(L); l(ji)>L. (18) 



In the thermodynamic limit, t(/i) diverges as r(/x) ~ e ( 16 / 3 ) vm*-m . Since \i ~ 1/T, r has a 
Vogel-Fulcher type divergence r(T) ~ e 



(16/3) 



VT-T, 



We have recently shown that the Vogel-Fulcher divergence of the dimer model can be 
obtained from an exact solution of the continuum version of the master equation if we 
assume that S(p) has a quadratic form jsi|. 



NUMERICAL ANALYSIS OF THE MASTER EQUATION 

We have carried out numerical diagonalizations of the rate matrix for LxL systems in 
order to verify some of the assumptions that have been made in the scaling analysis and 
in the calculation of the correlation function. These computations also provide us with 
information about the finite-size effects on the critical dynamics of the dimer model. 

The probability distribution, P(p,t), can be written in terms of the eigenvalues, A«, and 



eigenfunctions, ipi{p) of the rate matrix 23]: 



P(P,t) = E^PK^ (19) 

The eigenvalues of the rate matrix are non-negative and the equilibrium distribution 
is given by the zero-eigenvalue function, tpi (Ai = 0). The smallest, non-zero eigenvalue 
characterizes the state with the longest relaxation time. All correlation functions can be 
expressed in terms of the eigenvalue spectrum and, in particular, the equilibrium, tilt-tilt 
autocorrelation function can be written as: 

c(t) = Y. i ^Y l ^~ {Fp+F ^ip)^W (2°) 

Comparing to Eq. |H] it follows that the approximate form is obtained in the limit of delta- 
function localized eigenfunctions. We will show below that the eigenfunctions corresponding 
to non-zero eigenvalues of the dimer model are indeed well localized. 
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Results of the numerical diagonalization, using the exact entropy function, show that 
the longest time scale, r, of the dimer model increases in a non-arrhenius, Vogel-Fulcher 
fashion, as shown in Fig. 0] The scaling in this figure is what is expected from the scaling 



solution of the model with the quadratic entropy [31| and similar to the results obtained 
from the scaling arguments presented in this paper, however the length scale emerging is 
l(u) = [-p*zrj]} 2 n °t KlA = as wou ld be expected from the quadratic entropy results 

m 

The exact diagonalization results show, unambiguously, that the Vogel-Fulcher law char- 
acterizes the time scale divergence at the entropy vanishing transition in the dimer model. 



Eigenfunctions 

In the dimer model, the eigenfunction corresponding to the smallest non-zero eigenvalue 
is localized at the largest barrier, i.e, the largest value of p. Higher eigenfunction move to 
smaller barriers but are still localized. The expression we used for C(t) is exact for delta- 
function localization of eigenfunctions and the numerical results justify this assumption, a 
posteriori. 



Sensitivity of dynamics to barrier size 

The sensitivity of the relaxation times to the barriers heights has been 
investigated by using the quadratic entropy model and writing B(p) = 
e^(S(p-l/L)-S(p) + (E(p-l/L)-E(p)))/2] and varying c between and 1. 
For c = 0, we recover the usual Langevin dynamics and the time scales should increase as a 
power law and for c = 1 we have barriers corresponding to the loop dynamics. The results 
plotted in Fig |U] demonstrate that, for c = 0, r ~ (// — p)^ 1 which is consistent with a 
dynamical exponent z = 2 and a correlation length exponent of v = 1/2; the exponents 
expected from a Langevin description of a mean-field model. It is also clearly seen from this 
figure that even for c = 0.25, the timescale increases more rapidly than a power law. An 
analysis of the continuum limit of the dimer model dynamics shows that, for any non-zero 
value of c, r ~ e A ^ 1 ^ where A is a constant jy], 32]. These results taken all together 



imply that there is a whole class of systems, where dynamical constraints may lead to 



12 



0.1 



0.01 



0.001 =- 



0.000J 



i — i i i — i — i — r~ rrnn-| 1 — i — r~ i-rrrrj 1 — i — rrm 



♦ 



• L=4 
■ L=8 

♦ L=12 
A L=16 



01 



J I I I I I I I I I I 

0.1 1 



10 



100 



l([L)/L 



FIG. 4: Scaling of r(/x, L) in the dimer model. The figure shows that a scaling form can be 
constructed in terms of the lengthscale l{n) = (p* / (fx* — n)) 2 and the particular form demonstrates 
the Vogel-Fulcher scaling r ~ e^^ 2 . 

non-zero values of c, which belong to different universality classes of dynamical critical 
phenomena. These are characterized by a Vogel-Fulcher rather than a power-law divergence 
of relaxation time-scales, with c being an indicator of fragility [11]. In real systems, such 
as supercooled liquids, one expects that there is a large but finite energy scale at which 
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FIG. 5: Plots of the eighth and twenty-fifth eigenfunctions of the dimer model for L=16 and fi ~ /x*. 
These two eigenfunctions were chosen to illustrate the localization of the eigenfunctions and the 
shift towards p = with increasing spectral index 

the hard constraints are violated. This energy scale then leads to a long-time cutoff of the 
Vogel-Fulcher behavior. This time scale, may however, be well beyond any experimentally 
measurable time scales. 
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FIG. 6: Timescale r in the dimer model for different values of the barrier strength, c, plotted as a 
function of l(p) = f 

ADAM-GIBBS SCENARIO 

The entropy of the dimer model, S con f, which corresponds to S(p) evaluated at the 
equilibrium value of the order parameter p, goes to zero at the transition. Furthermore, our 
results clearly show that the longest time scale diverges in a Vogel-Fulcher manner. The 
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Adam-Gibbs relation, r = r e A ^ Sconf ^ T \ however, does not capture the physics since in the 
dimer model, S con f jumps from a finite value at p = to zero, as p changes discontinuously 
to 1. Thus, in this model, the exponential divergence of the relaxation time-scale at the 
transition is not accompanied by a continuous vanishing of S con f. The analysis presented in 
this paper clearly demonstrates that the Vogel-Fulcher divergence is rooted in the constrains 
which lead to loop dynamics. This type of nonlocal dynamics leads to the unusual transition 
rates with energy-lowering transitions being determined by changes in entropy and, therefore, 
to an exponential decrease of the number of energy-lowering trajectories as one approaches 
the zero-entropy state. 

The configurational entropy of supercooled liquids has been interpreted as the inherent 
structure entropy, i.e., the number of valleys at the temperature of interest. The observation 
that the Adam-Gibbs scenario describes much of the phenomenology of supercooled liquids 
could imply that there is a phase transition in the inherent-structure space similar to the 
one discussed in this paper for dimers. Experiments and simulations have shown that a hall- 
mark of supercooled liquids approaching the glass transition is the appearance of dynamical 
heterogeneities P. I33I. l^. The loops in the dimer dynamics are analogs of these dynamical 
heterogeneities since they define the correlated moves allowed by the constraints. These 
heterogeneities are present as long as the constraints are not violated and are characterized 



23,123. The anal - 



by a size distribution which changes as the critical point is approached 
ogy between loops and dynamical heterogeneities suggest that the dynamics in the inherent 
structure space of supercooled liquids could be similar to the loop dynamics of dimer models. 
If this is the case then transition rates between inherent structures should exhibit features 
similar to the ones discussed in this paper. We are currently in the process of analyzing 
transition rates between inherent structures of Lennard- Jones glass formers in order to get a 
better understanding of the connection between dynamical heterogeneities and the effective 
dynamics in the inherent-structure space. 
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